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ABSTRACT 

We present DES13S2cmm, the first spectroscopically-confirmed superluminous supernova 
(SLSN) from the Dark Energy Survey (DES). We briefly discuss the data and search algorithm 
used to find this event in the first year of DES operations, and outline the spectroscopic data 
obtained from the European Southern Observatory (ESO) Very Large Telescope to confirm 
its redshift (z = 0.663 ± 0.001 based on the host-galaxy emission lines) and likely spectral 
type (type I). Using this redshift, we find M^ eafc = —21.05lgjg for the peak, rest-frame 
(7-band absolute magnitude, and find DES13S2cmm to be located in a faint, low metallicity 
(sub-solar), low stellar-mass host galaxy (log(M/M Q ) = 9.3 ± 0.3); consistent with what is 
seen for other SLSNe-I. We compare the bolometric light curve of DES13S2cmm to fourteen 
similarly well-observed SLSNe-I in the literature and find it possesses one of the slowest 
declining tails (beyond +30 days rest frame past peak), and is the faintest at peak. Moreover, 
we find the bolometric light curves of all SLSNe-I studied herein possess a dispersion of 
only 0.2-0.3 magnitudes between +25 and +30 days after peak (rest frame) depending on 
redshift range studied; this could be important for ‘standardising’ such supernovae, as is done 
with the more common type la. We fit the bolometric light curve of DES13S2cmm with two 
competing models for SLSNe-I - the radioactive decay of 56 Ni, and a magnetar - and find 
that while the magnetar is formally a better fit, neither model provides a compelling match 
to the data. Although we are unable to conclusively differentiate between these two physical 
models for this particular SLSN-I, further DES observations of more SLSNe-I should break 
this degeneracy, especially if the light curves of SLSNe-I can be observed beyond 100 days in 
the rest frame of the supernova. 
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1 INTRODUCTION 

The last five years have seen the emergence of a new class of ultra- 
bright stellar explosions: superluminous supernovae (SLSNe; for 
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a review see lGal-Yamll2012l) . some 50 times brighter than classi¬ 
cal supernova (SN) types. Data on these rare and extreme events 
are still sparse, with only ~50 SLSN detections reported in the lit¬ 
erature. These events have generally been poorly studied with in¬ 
complete imaging and spectroscopy, leaving many aspects of their 
observational characteristics, and their physical nature, unknown. 

Yet these SLSNe could play a key role in many diverse ar¬ 
eas of astrophysics: tracing the evolution of massive stars, driv¬ 
ing feedback in low mass galaxies at high redshift, and providing 
potential line- of-sigh t probes of interstellar medium to their high- 
redshift ho sts i Berger et alJl2012l) . SLSNe have also been detected 
to 2 ~ 4 JCooke et al. 2012l~ far beyond the reach of the current 
best cosmological probe, type la supernov ae (SNe la). If SLSNe 
can be standardised (e.g. Insemi&Smartt| 201 4). as is done with 
SNe la dTripdl 19981; IPhillinjl 1 9951: iRiess et alJll996h . then a new 
era of SN cosmology would be possible, with the potential to accu¬ 
rately map the exp ansion rate of the universe far into the epoch of 
deceleration foelubac et al ■I20I4I) . 

SLSNe have been divided into three possib le types: SLSN-II, 
SLSN-Iand SLSN-R(seeFig. 1 ol lGal-Yanj20ll) . SLSNe-II show 
signs of interactions with CSM v ia narrow hydrogen lines (e.g., 
lOfek et alj2007l ;lSmit h et alj2007t) . and thus may simply represent 
the bright end of a continuum of Type Iln SNe (although this is not 
well established). 

SLSN-I a re sp ectroscopically classified as hydrogen free 
l lOu i mb v etalj 200^ and are p ossibly related to Type Ic SNe at 
late times iPastorello et al.l201()i) . but normal methods of powering 
such SNe (e.g., the radioactive decay of 56 Ni or energy from the 
gravitational collapse of a massive star) do not appear to be able to 
simultaneously reproduce their extreme bri ghtne ss, slowly-rising 
light curves, and decay rates dlnserra et aI] l2Q13h . Many alterna¬ 
tive models have been proposed: the injections of energy into a SN 
ejecta via the spin down of a young mag netar dKasen & BildsteiJ 
l20ld : IWooslevl l2010l : llnserra et alj 120131) : interaction of the SN 
ejecta with a massive (3-5 Mp,) C/O-rich circumstellar material 
(CSM; e.g., iBlinnikov & Sorokina! |2010|) ; or collisions between 
high-velocit y shells g enerated from a pulsational pair-instability 
event dWooslev et al. 120071) . These explanations are still actively 
debated in the literature. 

SLSNe-R are rare and characterised by possessing ex¬ 
tremely long, slow-declining light curves (> 200 days in the rest 
frame). SLSNe-R originally appeared consistent with the de ath of 
>100Mq stars via the pair instability mechanism iGal-Yam et alj 
120091 ). However, new observations - and the lack of significant 
spectral differences between SLSN-I and SLSN-R - ha ve c hal- 
lenged the notion that these classes are truly distinct dNicholl et alj 
uOO]) and maybe better described together as Type Ic SLSNe 
I Inserra et alJIbool) . 

In summary, the origin of the power source for SLSNe (of all 
types) remains unclear, with the possibility that further sources of 
energy may be viable. A more detailed understanding requires an 
increase in both the quantity and quality of the data obtained on 
these events. However, finding more SLSNe (of any type) with ex¬ 
isting transient searches is challenging, especially in the local uni¬ 
verse wheretherate^are only ~ 10~ 4 that of t he core-collapse SN 
rate dOuimbv et all2013l ; McCrum et al. l20l4) thus making it hard 
to assemble large samples of events over reasonab le timescales, 
even from large surveys. For example, iMcCrum et alj d2014l) de¬ 
tected only 10 SLSN candidates over 0.3 ^ z ^ 1.4 within the 
first year of the PanSTARRS 1 Medium Deep Survey. 

Searches for SLSNe at higher redshifts may be more prof¬ 
itable as the rates of SLSNe appear to rise by a factor of ~10- 


15 at z > 1.5 dCooke et al.ll2012l) . perhaps tracking the increased 
cosmic star-formation rate and/or decreasing cosmic metallicity: 
SLSNe are preferentially found in faint, low-metallic ity galaxies 
jNeill et al.l201ll : IChen et al.l2013l : lLunnan et all2014l) . Moreover, 
the peak of a SLSN’s energy output is located in the UV region of 
the electromagnetic spectrum, which is redshifted to optical wave¬ 
lengths at high redshift. 

In this paper, we outline our first search for SLSNe in the 
Dark Energy Survey (DES). In Section [2] we discuss details of 
DES and our preliminary search, while in Section [3] we describe 
the first SLSN detected in DES and provide details of its proper¬ 
ties, including our spectroscopic confirmation. In Section [4] we 
compare our SLSN to other events in the literature and explore 
the possible power sources for our event. We conclude in Sec¬ 
tion [5] Throughout this paper, we assume a flat A-dominated cos¬ 
mology with Q.m = 0.28 and Ho = 70kms _1 Mpc - 1 , con sistent 
with recent cosmolo gical measurements (e.g. lAnderson et al J2014I: 
iBetoule et alj|2014h . 


2 THE DARK ENERGY SURVEY 


The Dark Energy Survey (DES; Elaueheill2005i ) is a new imag¬ 
ing survey of the southern sky focused on obtaining accurate con¬ 
straints on the equation-of-state of dark energy. O bservatio n s are 
carried out using the Dark Energy Came ra (DECam ; Flaugher et al. 


120121; Diehl & For the Dark Energy Survey Coll aborati on! 120121) 
on the 4-metre Blanco Telescope at the Cerro Tololo Inter- 
American Observatory (CTIO) in Chile. DECam was successfully 
installed and commissioned in 2012, and possesses a 3 deg 2 field- 
of-view with 62 fully depleted, red-sensitive CCDs. 


2.1 DES Supernova Survey 


DES conducted science verification (SV) observations from 
November 2012 to February 2013, and then began full science op¬ 
erations in mid-August 2013, with year one running until Febru¬ 
ary 2014. DES is performing two surveys in parallel: a wide-field, 
multi-colour ( grizY ) survey of 5000 deg 2 for the study of clus¬ 
ters of galaxies, weak lensing and large scale structure; and a 
deeper, cadenced multi-colour ( griz ) search for SNe. In this pa- 
per we focus on this DES SN Survey, as outlined in lBernstein et alj 
d2012l) . which surveys 30 deg 2 over 10 DECam fields (two ‘deep’, 
eight ‘shallow’) located in the XMM-LSS, ELAIS-S, CDFS and 
‘Stripe82’ regions. As of November 2014, the DES SN Survey has 
already discovered over a thousand SN la candidates, with more 
than 70 spectroscopic classifications to date. 

This rolling search for SNe la is also ideal for discovering 
other transient phenomena including SLSN e. Us ing the light curve 
of SNLS-06D4eu (a SLSN-I; iHowell et alj2013l) . we calculate that 
DES could detect such an event to z = 2.5 given the limiting mag¬ 
nitudes for the DES deep fields. For example, assuming a volumet¬ 
ric rate of 32 Gpc~ 3 yr _1 from lOuimbv et all d2013l) and correct¬ 
ing for time dilation, we estimate DES should find approximately 
3.5 SLSN-I events per DES season (five months) in the deep fields 
and many more in the shallow fields (although at a lower redshift). 
This basic prediction is uncertain because: i) We have not corrected 
for ‘edge effects’ that will reduce the number of SLSNe with sig¬ 
nificant light-curve coverage; ii) The volum etric r ates themselves 
may evolve with redsh ift (e.g., Cooke et~ahll20 1 2|) and have large 
uncertainties ( IMcCrum et al.l!20T4h : and iii) We have not included 
any corrections for survey completeness. 
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2.2 Selecting SLSNe candidates 


3.1 Light Curve 


We began our search for SLSNe using the data products from the 
DES SN Survey. All imaging data were de-trended and co-added 
using a standard photometric reduction pipeline at the National 
Center for Supercomputing Applica tions (NCSA) using the DES 
Data Management system fDESDM: lMohr et aPl2012l : lDesai et ail 
|20I3) . producing approximately 30 ‘search images’ per field (in all 
filters) over the duration of the five-month DES season. We perform 
difference imaging on each of these search images, using deeper 
template images for each field created from the co-addition of sev¬ 
eral epochs of data obtained during the SV period in late 2012. Be¬ 
fore differencing, the search and template images were convolved 
to the same point spread function (PSF). 

Objec ts were selected from the difference images using SEx- 
TRACTOR llBertin & Arnoutslll996l) v2.18.10, and previously un¬ 
known transient candidates were identified and examined using vi¬ 
sual inspection (or ‘scanning’) by DES team members. Over the 
course of the first year of DES approximately 25, 000 new tran¬ 
sient candidates were selected in this manner, including many data 
artefacts in addition to legitimate SN-like objects. A machine learn¬ 
ing algorithm (Goldstein et al., in preparation) was used to improve 
our efficiency of selecting real transients. Transient candidates were 
photometrically cl assifi ed using the Photometric SN IDentification 
software (PSNID; ISako et af3l2008l 1201 ll) to determine the likely 
SN type based on fits to a series of SN templates. These SN can¬ 
didates were stored in a database and prioritized for their spectro¬ 
scopic follow-up. 

During the first season of DES, PSNID only included tem¬ 
plates for the normal SN types Ib/c, II and la, so a fully automated 
classification of SLSN candidates was not possible. Therefore, for 
this first season of DES, we searched for candidate SLSNe using 
the following criteria: i) At least one month of multi-colour data, 
i.e., typically five to six detections ( S/N > 3.5) in each of griz\ 
ii) A low PSNID fit probability to any of the standard SN sub¬ 
classes; iii) Located greater than one pixel from the centroid of the 
host galaxy (to eliminate AGN); and iv) Peak observed brightness 
no fainter than one magnitude below that of its host galaxy. These 
broad criteria will be satisfied by most SLSNe, while also helping 
to eliminate many of the possible contaminating sources. 

Using this methodology, we selected 10 candidate SLSNe 
over the course of the first season. We secured a useable spectrum 
for one of our SLSN candidates, DES 13S2cmm, which is discussed 
in detail in Section[3] while none of the other candidates were spec¬ 
troscopically confirmed. We are attempting to obtain spectroscopic 
redshifts from the host galaxies of these other SLSN candidates, 
though this will be challenging in the several cases where no host 
is detected in our template images. 


3 DES13S2CMM 

In Fig. |T] we show DES riz imaging data at the position of 
DES13S2cmm: pre-explosion, post-explosion, and their difference 
image. DES13S2cmm is located at RA(2000) = 02 h 42 m 3Z82, 
Dec. (2000) = —01°21 , 30 , /l, and was first detected on 2013 Au¬ 
gust 27 in the DES SN ‘S2’ field (a ‘shallow’ field) located in the 
Stripe 82 region. DES transient names are formatted following the 
convention DESYYFFaaaa, where YY are the last two digits of the 
year in which the observing season began (13), FF is the two char¬ 
acter field name (S2), and the final characters (all letters, maximum 
of 4) provide a running candidate identification, unique within an 
observing season, as is traditional in SN astronomy (cmm). 


In Fig.[2]we present the multi-colour light curve for DES 13S2cmm 
constructed from the first year of DES observations. We have not 
included upper limits from the earlier, pre-explosion epochs avail¬ 
able from the DES SV period in 2012-2013, as these data were 
taken hundreds of days prior to the data shown in Fig. [2] and are 
already used to create the reference template image for this field 
used in the difference imaging detection of the_SNe. 

We use the photometric pipeline from ISullivan et al .1 d201 lh 
(also used in lMaeuire et alj|2()l2l : l0fek et all 2013 ) to reduce the 
DES photometric data for DES 13S2cmm. This pipeline is based on 
difference imaging, subtracting a deep, good seeing, pre-explosion 
reference image from each frame the SN is present in. The pho¬ 
tometry is then measured from the differenced image using a PSF- 
fitting method, with the PSF determined from nearby field stars in 
the unsubtracted image. This average PSF is then fit at the position 
of the SN event, weighting each pixel according to Poisson statis¬ 
tics, yielding a flux, and flux error, for the SN. 

We calibrate the flux measurements of DES13S2cmm to a set 
of tertiary standard stars produced by the DES collaboration (Wyatt 
et al., in preparation ), intended to be in the AB photometric system 
( joke & Gunnl 19831) . The photometry shown in Fig.[2]has also been 
checked against the standard DES detection photometry, used to 
find and classify all candidates for follow-up, and found to be in 
good agreement. The light-curve data are presented in Tablc|T| 

We correct for Galac tic extinction using the maps of 
ISchlaflv & Finkbeineil d201 ll) . which estimate E(B — V) = 0.028 
at the location of DES13S2cmm, leading to extinction values of 
0.123, 0.070, 0.051 and 0.039 magnitudes in the DES griz filters 
respectively. We do not include these corrections in the values re¬ 
ported in Table |T] but do correct for extinction in all figures and 
analysis herein, including Figure[2] 


3.2 Spectroscopy 

On 2013 November 12, we requested Director’s Discretionary 
Time at the European Southern Observatory (ESO) Very Large 
Telescope (VLT) to observe DES13S2cmm, and were awarded 
target-of-opportunity time within days. The closeness of the ob¬ 
ject to the moon and instrument scheduling issues conspired 
to delay observations until 2013 November 21, at which point 
DES13S2cmm was approximately 30 days after peak brightness 
(rest frame). A spectrum was obtained with the FOcal Reducer and 
low dispersion Spectrograph (FORS2 fAppenzeller et alJll998l) us¬ 
ing the GRIS_300I+11 grism, the OG590 order blocker, and a 1" 
slit, with an exposure time of 3600s (3x 1200s). This configuration 
provided an effective wavelength coverage of 5950 A to 9400 A. 

The reduction of the FORS2 data followed standard proce¬ 
dures using the Image Reduction and Analysis Faci lit\[] environ¬ 
ment (v2.16), using the pipeline described in lEllis et alJ J20081 ). 
This pipeline includes an optim al two-dimensio nal (2D) sky sub¬ 
traction technique as outlined in lKelsoiil 1200 3l) . subtracting a 2D 
sky frame constructed from a sub-pixel sampling of the background 
spectrum and a knowledge of the wavelength distortions deter¬ 
mined from 2D arc comparison frames. The extracted ID spectra 
were then scaled to the same flux level, the host-galaxy emission 


1 The Image Reduction and Analysis Facility (IRAF) is distributed by the 
National Optical Astronomy Observatories, which are operated by the As¬ 
sociation of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Figure 1. Colour images (riz) of the field surrounding DES13S2cmm. Left: Deep template image created from the co-addition of several epochs of data obtained 
during the DES Science Verification period in late 2012; the likely host galaxy is indicated by the arrow. Centre: The search image taken close to maximum 
light (28-Sept-2013) with the position of DES13S2cmm indicated. Right: The difference of the two previous images, clearly identifying DES13S2cmm at 
center. This image demonstrates the quality of difference images used for discovery and monitoring during the DES observing season. 


1-Sep-2013 1 -Oct-2013 1-Nov-2013 1-Dec-2013 1-Jan-2014 1-Feb-2014 



Figure 2. The observed AB magnitude light curve of DES13S2cmm in the four DES SN search filters ( griz ) as a function of the observed phase. ‘S’ denotes 
the epoch when the ESO VLT spectrum was obtained. The error-bars represent 1-cr uncertainties, and the arrow symbols are 3-cr upper limits. The phase 
relative to peak was defined using the observed-frame r-band. We have artificially offset (in magnitudes) the r, i and 2 band data for clarity. 



Figure 3. The reduced VLT FORS2 two-dimensional spectrum for DES13S2cmm. The key nebular emission lines used to derive the redshift of the host 
galaxy, and therefore supernova, are labelled. 
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Table 1. Light curve data for DES13S2cmm used in Fig. [2] We provide the observed calendar date, modified Julian Day (MJD), observed phase relative to 
peak MJD (defined relative to peak flux in the r-band), and fluxes (not corrected for Galactic extinction) with 1-cr uncertainties in the DES griz passbands. 
These fluxes can be converted to AB magnitudes using mag AB = —2.51og 10 (/) + 31. 


Date 

MJD 

Phase 

(days) 

fa 

fr 

fi 

fz 

30-AUG-2013 

56534.3 

-25 

1138+77 

1705+94 

1730+181 

1913+231 

3-SEP-2013 

56538.3 

-21 

1112+103 

2063+129 

2234+186 

1845+288 

8-SEP-2013 

56543.3 

-16 

1369+104 

2348+148 

2285+198 

2285+277 

12-SEP-2013 

56547.2 

-12 

1373+130 

2329+167 

2433+207 

2376+284 

15-SEP-2013 

56550.2 

-9 

± 

2691+413 

2818+380 

2193+389 

24-SEP-2013 

56559.2 

0 

2022+388 

3727+363 

3183+346 

3548+295 

28-SEP-2013 

56563.2 

4 

1668+97 

3688+133 

3400+187 

3488+241 

2-OCT-2013 

56567.2 

8 

1805+141 

3666+166 

3742+209 

3169+261 

10-OCT-2013 

56575.2 

16 

1370+108 

3285+137 

3636+192 

3182+230 

14-OCT-2013 

56579.1 

20 

1260+243 

3182+223 

3482+259 

2663+295 

25-OCT-2013 

56590.3 

31 

886+346 

1782+319 

2924+277 

3163+352 

29-OCT-2013 

56594.1 

35 

667+146 

1782+187 

2398+267 

2987+366 

6-NOV-2013 

56602.1 

43 

344+103 

1698+132 

2326+185 

2433+257 

10-NOV-2013 

56606.1 

47 

803+426 

1338+418 

2633+510 

1817+634 

13-NOV-2013 

56609.1 

50 

-312+276 

970+252 

2672+281 

2726+426 

20-NOV-2013 

56616.1 

57 

347+195 

1855+239 

1848+307 

1739+331 

I-DEC-2013 

56627.0 

68 

184+107 

1306+129 

2032+176 

2322+216 

9-DEC-2013 

56635.0 

76 

259+208 

1322+178 

1766+225 

2318+304 

19-DEC-2013 

56645.1 

86 

302+226 

1305+211 

1697+230 

1808+248 

23-DEC-2013 

56649.1 

90 

261 + 136 

956+157 

1667+248 

2038+276 

27-DEC-2013 

56653.1 

94 

41 + 156 

1295+174 

1613+201 

1930+248 

3-JAN-2014 

56660.0 

101 

252+119 

1000+132 

1736+196 

1993+261 

10-JAN-2014 

56667.1 

108 

310+343 

1144+364 

1609+408 

1586+387 

18-JAN-2014 

56675.0 

116 

-75+236 

717+220 

1290+261 

1863+297 

25-JAN-2014 

56682.0 

123 

306+158 

992+188 

1382+215 

1254+270 

30-JAN-2014 

56687.0 

128 

345+188 

523+179 

1518+217 

1404+261 

6-FEB-2014 

56694.0 

135 

-346+359 

1459+376 

1602+465 

2443+461 


lines interpolated over, and the spectra combined using a weighted 
mean and the uncertainties in the extracted spectra. 

In Fig [3] we present our reduced 2D FORS2 spectrum, 
where wavelength increases to the right along the horizontal axis 
and the vertical axis is distance along the slit. We highlight 
the obvious narrow nebular emission lines in the host galaxy 
of DES13S2cmm, specifically [Oil] A3727A, H/3 A4861A, and 
[O III] A4959, 5007 A. We assume these lines originate from the 
underlying host galaxy as their spatial extent along the vertical axis 
is greater than the observed trace of the main supernova spectrum. 
This observation does not exclude the possibility that such emis¬ 
sion lines come from the supernova, but it is difficult to address 
this issue with our low resolution, and low signal-to-noise, data. We 
note the observed emission lines are not significantly broadened or 
asymmetrical, as witnessed in type Iln supernovae. 

We see no evidence of any further emission lines, from the SN 
and/or host galaxy. Therefore, it is unlikely DES 13S2cmm is a type 
II supernova (normal or superluminous). 

Using these nebular emission lines, we determine a redshift of 
z = 0.663 ± 0.001 for the host galaxy of DES13S2cmm, where 
the uncertainty is derived from the wavelength dispersion between 
individual emission lines. This measurement is consistent with the 
redshift of the SN based on identifications of the broad absorption 
features seen in the spectrum. As the redshift derived from the host- 
galaxy spectrum is significantly more precise, henceforth we adopt 
this value as the redshift to DES13S2cmm. 

The spectral classifi cation for DES13S2cmm is based on the 
SUPERFIT program dHowell et al .120051) . which compares observed 


spectra of SNe to a library of template SN spectra (via \ 2 rnin- 
imisation) while accounting for the possibility of host-galaxy con¬ 
tamination in the data. We added as template spectra of SLSNe- 
I (PTF09atu, PTF09cnd, PTFlOnmn) and SLSN-R (SN2007bi) 
to SUPERFIT to facilitate our classification. We found that the 
highest-ranked fits to DES13S2cmm were these SLSNe templates, 
preferred over templates of normal SN types (e.g., SN la). The 
SLSNe-I templates were slightly better fits to DES13S2cmm than 
SN 2007bi (SLSN-R), although this may be due to not having a 
template for SN 2007bi at a similar phase to our spectrum for 
DES13S2cmm. 

We show our spectrum of DES13S2cmm in Fig. [4] alongside 
spectra from other well-st udied SLSNe in the l i terature at a similar 
phase in their light c urves dOuimbv et alj|201~ll : llnserra et alj|2013l : 
iGal-Yam et al 1200(d) . There is a deficit of spectra for SLSNe in the 
literature at these late phases for a robust comparison. However, the 
broad resemblance between the spectra of PTFlOhgi, SN2011kf, 
SN 201 Ike and DES13S2cmm, all taken at approximately +30 
days past peak, suggests DES13S2cmm is a similar object to these 
other SLSNe-I. The spectrum of DES13S2cm m is also similar to 
the spectrum of SN2007bi (a SLSN-R in iGal-Yam et~al] I2009I) . 
although the differences between Type I and Type R SLSNe re¬ 
mains unclear, and either classification would be interesting given 
the sparseness of the data. We also show in Fig. [4] the spectrum 
of SN 2007if, a ‘super-Chandrasekhar" SN la dScalzo et al.ll2Qlol) . 
an over-luminous ( My = —20.4) and likely thermonuclear SN 
located in a faint host galaxy. This comparison clearly shows that 
DES13S2cmm is unlikely to be a super-Chandrasekhar SN la. 
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Rest wavelength (Angstroms) 

Figure 4. The VLT FORS2 spect rum of DES13S2c mm (second from top) com pared to other SL SNe of similar pha se in their light curves. From top to 
bottom, we also show SN2011kf dlnserra et alj|2013h , PTFlOhgi iOuimby et alj|20Tll) . SN 201 Ike ilnserra et alj|2013h and SN 2007bi(classified as Type R 
in Gal- Yam et alj|2009h ; as discussed in Section ns SUPERFIT identifies these SLSNe as some of the best-fitting spectral templates to DES13S2cmm. In 
each case, the light grey shows the raw data and the solid line is a Savitsky-Golay smoothed version of the spectrum. The spectra are labelled with the SN 
name, SLSN type (except DES13S2cmm), re dshift, and the phase when the spectrum was taken. We show the spectrum of the ‘super-Chandrasekhar’ SN la 
SN 2007if as discussed in Scalzo et al.1 l20ld) . Note that the spectrum for 2007bi, while at a later phase than the other data, is the earliest spectral observation 
of this SN. 


3.3 Host-galaxy properties 

The host galaxy of DES13S2cmm is detected in our reference im¬ 
ages that contain no SN light. Using SEXTRACTOR in dual-image 
mode and the i-band image as the detection image, we measure 
host-galaxy magnitudes of g = 24.24 ± 0.13, r = 23.65 ± 0.10, 
i = 23.31 ± 0.10 and z = 23.26 dh 0.17 (AB MAG-AUTO magni¬ 
tudes). The large photometric errors are due to the relatively shal¬ 
low depth, and poor seeing, of the template image from the DES S V 
data - we cannot use the DES first-year data for the host photometry 
as all epochs of this data contain some light from DES13S2cmm. 

The estimated photometric redshift of the host galaxy at the 
time of discovery was %,hoto = 0.86 ± 0.13 (68 percentile er¬ 
ror). This photometric redshift was derived by the DESDM neu¬ 
ral network photo-z module code, with uncertainties estimated 
from the nearest neighbour method, using data taken during SV 
dSanchez et al.ll20l4l) . Photometric redshifts (when available) are 
used when selecting SLSN candidates for spectroscopic follow¬ 
up, and for DES13S2cmm this redshift implied the event was 


brighter than M* eak = -21 (typical for SLSNe; Ic.al-Yandl2()12l) . 
Since the start of the first DES observing season the photomet¬ 
ric catalog derived from SV data has improved, and presently the 
DESDM photo-z for the host of DES13S2cmm has improved to 
2 P hoto = 0.71 ± 0.06, which is in better agreement with our mea¬ 
sured spectroscopic redshift (z = 0.663 ± 0.001). 

Using the DES host-galaxy photometry and the spectroscopic 
redshift from VLT, we estimate that the stellar mass of the host 
galaxy is log (M/Mq) = 9.3 ± 0 .3 using stellar popu lation 
models from PEGASE.2 dFioc & Rocca-Volmerangel f 1 997t) . and 
log( M/Mp)) - - 9.0 ± 0.3 using the stellar population templates 
from iMarastonl d2005lf (both corrected for Milky Way extinction 
and fixing the redshift of the templates to z — 0.663). Such a low 
stellar- mass host galaxy is consistent with the findings for o ther 
SLSNe dNeill et all201ll : IChen et all2013l:lLunnan et all 20 ill . 

We estimate the host-galaxy metallicity from our VLT spec¬ 
trum using the double-valued metallicity indicator R 23 , defined 
as ([Oil] A3727 + [OIII] AA4959, 5007)/T//3. Fluxes and un- 
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certainties are derived individually for each emission line from 
fits of a gaussian plus a first-order polynomial, with the assump¬ 
tion that there is no contamination of the galaxy emission lines 
from the SN. We compute m etall icities usin g the c alibrations of 
iKobulnickv & Kewlevl 1 20041) and McGauglJ 1 199ll) . and use the 
formulae derived in iKewlev & Ellison| j 2008 ) to convert these 
metallicities into the calibration of Kewlev & Dopital ( 120021 ) - a 
step that allows for direct comparison between different methods. 
We note that since Ha and Nil A6584 are redshifted to wave¬ 
lengths beyond our spectral coverage, we cannot determine the 
branch of the R 23 function that the metallicity lies on. However, 
the derived value for R 23 is close to its theoretical maximum, which 
minimizes the difference in metallicity estimates between the two 
branches. 

Assuming the lower (upper) branch , we find a metallic¬ 
ity (12 + logfO/ Hl) of 8.30 (8.38) from iMcGaued dl99ll) and 
8.30 (8.42) from IKobulnickv & Kewlevl d2004l) . which is consis¬ 
tent with the median metallicity of 8.35 found by iLunnan et alj 
d2014l) for a sample of 31 SLSNe host galaxies. The total un¬ 
certainty is dominated by the calibration (~ 0.15 dex), while 
the sta tistic al un certainty and the scatter induced by the conver¬ 
sion to [Kewlev & Dopita i d20021) are comparatively negligible. We 
thus find a host-g alaxy abundance ratio that is distinctly sub¬ 
solar (8.69; lAsplund et al.1 |2009I ). though not as extremely low 
as has been see n for other S LSNe, such as SN2010gx (12 + 
loglO/Hl=7.5: IChen et all2013l) . 


4 RESULTS 

4.1 Bolometric light-curve of DES13S2cmm 

We show in Fig. [5] the bolometric light curve of DES13S2cmm. 
This single light curve was constructed by fitting a single black- 
body spectrum to the DES multi-colour data (requiring measure¬ 
ments in a minimum of three photometric bandpasses for a fit), 
on each epoch in the light curve (Table 1). In detail, we redshift 
black-body spectra to the observer frame and integrate these spec¬ 
tra through the response functions of the DES filters. We then deter¬ 
mine the best-fitting parameters for a black body to the extinction- 
corrected photometry and integrate this spectrum in the rest frame 
to obtain the bolometric luminosity per epoch. This procedure ben¬ 
efits from the homogeneity of the DES data as we possess approx¬ 
imately equally-spaced epochs across the whole light curve - typi¬ 
cally in four DES passbands - due to our rolling search. The errors 
on the bolometric luminosities were calculated from the fitting un¬ 
certainties on the black-body spectra. 

We visually checked our best-fit blackbody curves for each 
epoch of the DES13S2cmm light curve against the observed griz 
photometry and found them to be reasonable within the errors. The 
largest discrepancies were seen in the g-band where some epochs 
were below the fitted curves, possibly due to UV absorption in 
the underlying spectral energy distributions relative to a blackbody 
curve. Inserra et al. (2013) compared the predicted UV flux (from 
a blackbody curve) for two SLSNe with SWIFT UV observations, 
and found the fluxes were consistent within the error. Therefore 
we make no additional correction to the derived bolometric light 
curve, but note that quantifiable systematic errors in the bolometric 
luminosity require spectral templates, which in turn requires more 
spectral data of SLSNe than is currently publicly available. 

Using the bolometric light curve of DES13S2cmm, we can 
determine the peak luminosity of the event and compare it to 


the definition of a SLSNe in iGal-Yaml I2OI2I) . This was achieved 
by integrating the best-fitting black-body spectra through a stan¬ 
dard (7-band filter response function (on the Vega system) to ob¬ 
tain the rest-frame absolute magnitude in (7-band (Mu)- We find 
the 2013 September 24 epoch (t = 0 relative phase in Table 1) 
gives the largest (peak) absolute ?7-band magnitude of M^ eak = 
—21.05lQ'Jg, which is consistent with the SLSNe threshold of 
M peak < 21 [ n | Gal _ Y aml i2012l) . 

We check this result using a model-independent estimate 
based on the observed r-band peak magnitude. At the redshift of 
DES13S2cmm, the observer-frame r band maps into the rest frame 
as a synthetic filter with an effective wavelength of « 3800A, 
slightly redward of the standard (7-band. Defining a synthetic filter 
as such means the k-correction is independent of the SN spectral 
energy distribution, and since we are in the AB system the correc¬ 
tion is simply 2.5 log(l + z). We apply this k-correction to the ob¬ 
served r -band peak magnitude, as well as corrections for Galactic 
extinction and distance modulus, to yield an estimated peak abso¬ 
lute magnitude of —20.47 (AB), or —20.74 (Vega) in our synthetic 
filter (where we have computed the AB to Vega conversion for 
our filter). While this is fainter than our blackbody-based estimate, 
Vega magnitudes are sensitive to the exact shape of the filter in this 
region due to the jump in Vega from 3700-3900A; thus even the 
small difference in effective wavelength between (7 band and our 
synthetic filter can result in significant offsets. The absolute magni¬ 
tude in our synthetic filter from the best-fitting blackbody method 
yields M synt hetic = —20.59, proving that our model-independent 
estimate is consistent with our results from the blackbody fitting. 

We note that the peak magnitude (Vega) for DES13S2cmm 
classifies it as a SLSN, provided the magnitude is intrinsic to the 
SN and not enhanced by other effects. Visual inspection of the im¬ 
ages in Fig. Q] (including the surrounding areas) shows little evi¬ 
dence for strong gravitational lensing which could affect the ob¬ 
served brightness of this SN (we note that the host galaxy would 
have to be lensed as well). We cannot conclusively rule out the pos¬ 
sibility of strong lensing, as in the case of PS 1 -1 Oafx iOuimbv et alj 
I2014I) . but as discussed above the spectrum and extended light 
curve of DES13S2cmm are inconsistent with a normal SN, while 
the probability of such a strong lensing event in our DES data is 
less than the probabilit y of di scover in g a S LSN, based on the lens¬ 
ing statistics given in Quimby et alj 1 120141) . and the lower redshift 
of DES13S2cmm and smaller search area of DES. Moreover, the 
observed r — i peak colour of DES13S2cmm is consistent with 
the expected colours of other unlen sed S Ne (including SLSNe) 
shown in Figure 4 of I Quimby et alj d201-4l) . i.e., DES13S2cmm is 
located below the thick black line in this plot. Therefore, we do 
not consider gravitational lensing further in this paper, but higher- 
resolution imaging data should be obtained for DES13S2cmm, af¬ 
ter the SN event has vanished, to provide better co nstra ints on pos¬ 
sible strong lensing as discussed in lOuimbv et alj d2014l) . 


4.2 Comparison of bolometric light-curves 

For comparison with DES13S2cmm, we also show in Fig. [5] 
the bolometric light cur ve s for fourteen SLSNe-I (or Ic, as 
discussed in Inserra etaL 2013) in the literatur e: PSl-lOky. 
PSl- lOawh dChomiuk et alj 1201 ll) and PSl-lObzj (Lunnanetal. 
120131) : SN201Qgx dPastorello et al1l2010l) : P TF lOhgi, SN201 Ike , 
PTFllrks, SN2011kf and SN20 12il dlnserra et alj 120131) : 
SNLS 06D4eu and SNLS07D2bv dHowell et al.l 120131); an d 
LSQ 12dlf, SSS 120810 and SN 2013dg dNicholl et alj 120141) . 
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Figure 5. The bolometric light curve of DES13S2cmm (red circles) compared to other SLSNe in the literature (see text for explanation). These bolometric light 
curves were constructed by fitting a single black-body spectrum to the multicolour data available per epoch (requiring at least three photometric measurements 
per epoch), and then by integrating the black-body spectral energy distribution to produce the bolometric luminosity per rest-frame epoch (described in detail 
in Section EES. We show statistical errors on the blackbody fits to each epoch of these light curves, but in most cases this error is smaller than the plotting 
symbol. We stress that these errors do not account for systematic uncertainties involved in the bolometric calculation (e.g., assuming a blackbody for the 
spectral energy distribution of SLSNe, especially for the UV part of the rest-frame spectral energy distribution). The phase relative to peak (y-axis) for the 
literature SLSNe is taken directly from the published individual studies and not recalculated herein. The dashed lines shown in both panels represent the 
expected *®Co decay rate at late times. In the lower panel, we show the bolometric light curve expected for a typical Type la SNe (solid line) to again illustrate 
the extraordinary brightness, and extended duration, of these SLSNe 


We exclude from this figure, and our analysis, other SLSNe-I 
that do not possess data for at least three passbands for each 
epoch (defined as being taken within a 24 hour window) over a 
majority of their light curves. This criterion, so defined to ensure 
well-constrained bolometric light-curve fits, excludes namely 


PTF09atu, PTF09cnd, and PTF09cwl d Ouimbv et al J 2011): 
SCP06F6 dBarbarv et aid [200(d) : and SN2006oz i Leloudas et aid 

SqU). 


The full bolometric light curves for these SLSNe-I were 
calculated using the same methodology as described above for 
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DES13S2cmm (i.e. using the published photometric data and re¬ 
quiring three or more passbands per epoch). As a check, our 
bolometric light curves for PTF lOhgi, SN 201 Ike, PTFllrks, 
SN2011Jcf and SN2012il are similar to those published in 
llnserra et all d2013l) . with any offsets due to differences in the inte¬ 
gration o f the blac kbody curve; we integrate over all wavelengths, 
while Inserra et alj J 2 OI 3 I) integrate over the wavelength range of 
their observed optical passbands (we are able to reproduce their 
luminosities if we follow their procedure). Our bolometric light- 
curves for PSl-lOky. PSl-lOawh, and PSl-lObzj are also similar 
to those published in the literature. 

In Fig. [5] the phase relative to peak for DES13S2cmm was de¬ 
fined as the peak of the observer-frame r-band, which occurred on 
MJD 56563.2, while for other SLSNe-I we simply used the phase 
relative to peak as reported by the individual literature studies. In 
the case of SNLS 06D4eu and SNLS 07D2bv, iHowell et all d20 1 3l) 
reported a phase based on their calculation of the bolometric light 
curves, which we use herein, but note there may be greater uncer¬ 
tainty in the definition of the time of peak for these light curves. 

The top panel of Fig. [5] shows that DES13S2cmm is one of 
the faintest SLSNe-I around peak compared to the other SLSNe- 
I in the literature, and possesses one of the slowest declining 
tails (beyond approximately +30 rest-frame) of the SLSNe studied 
herein. In the bottom panel of Fig. [5] we compare DES13S2cmm 
to possible Type R SLS Ne (to be consistent with the classifica¬ 
tion scheme of lGal- Yam! 20121) and find the l ate-time li ght curve of 
DES13S2cmm is similar to both PSl-llap dMcCrum et al.ll2()l4h 
and SN2007bi dGal-Yam et alj 120091) . highlighting the possible 
overlap between th ese t w o S LSNe types, as discussed in Sec- 
tion l3.2l and Inserr a et aid d2013l) . 

Fig. [5] also demonstrates that the fifteen SLSNe-I light curves 
(top panel) have similar luminosities at approximately +25 days 
(rest frame) past peak, close to the inflection point in the light 
curves when their extended tails begin to appear. We calculate the 
dispersion, as a function of relative phase, for the SLSNe-I light 
curves shown in Fig. [5] linearly interpolating between large gaps 
in each light curve to provide an evenly sampled set of data. We 
find a minimum dispersion of cr(logioLboi) = 0.11 - equivalent to 
0.29 magnitudes - around +30 days past peak (rest frame), derived 
from the nine SLSNe in Fig.[5]constraining this phase. This disper¬ 
sion in the light-curves can be decreased to cr(logioLboi) = 0.083, 
or 0.20 magnitudes, at approximately +25 days past peak (rest 
frame) if we exclude the four highest redshift SLSNe-I above a ~ 1 
(SNLS06D4eu, SNLS07D2bv, PSl-lOky. PSl-lOawh). The bolo¬ 


metric luminosities of these SLSNe-I are more difficult to esti¬ 
mate because of the uncertainties associated with modelling the 
UV part of the spectrum by assuming a simply blackbody curve. 
There could al so b e evolutio n in the population with redshift. As 
commented in lOuimbv et al.l ( |201 lb . such light curve characteris¬ 
tics could point towards a possible ‘standardisation’ of SLSNe-I, 
as with SNe la, leading to their u se as hig h-redshift ‘standard can¬ 
dles’ for cosm ological studies fang et all20 f+l ). 

Recently, llnserra & Smartti 1 2014i) has proposed a ‘standard¬ 
isation’ of the peak magnitude of 16 SLSNe using a stretch- 
luminosity correction similar to the Amis relationship of IPhillipsI 
d 19931) for SNe la. They find that SLSNe that are brighter at peak 
have slower declining light curve over the first 30 days after peak 
(rest frame) and, correcting for this correlation, the dispersion of 
the peak magnitude of SLSNe can be reduced to a ~ 0.2. This is 
similar in size to the dispersion between +25 to +30 days past peak 
seen in Fig. [5] and is only a factor of two worse than the dispersion 
obtained for standardised SNe la. 


Both Fig. [5l and llnserra & Smartti J 20141) suggest some stan¬ 
dardisation of SLSNe light curves is possible but further obser¬ 
vations and analyses are required to find the best approach. For 
example, the light curve of D ES13S2cmm does not a ppear to fol¬ 
low the correlation found by llnserra & Smartl) J201411 as it has a 
faint peak magnitude, relative to other SLSNe, but is also slowly 
declining towards +30 days after peak. A direct comparison to 
llnserra & Smartti 120141) is difficult as our study uses different sam¬ 
ple selection criteria (only 9 of their 16 SLSNe are included here 
becau se of our requirement t o have three passbands per epoch), 
while llnserra & Smartll ( l2014h estimate a common synthetic peak 
magnitudes for all their SLSNe using their own time-evolving spec¬ 
tral template. A more detailed comparison will be possible with 
further DES SLSNe. 


4.3 Power source of DES13S2cmm 

The details of the power source of SLSN-I remain unclear and are 
much debated (see Section |T}. Here we investigate the two popu¬ 
lar explanations in the literature for SLSNe-I: radioactive decay of 
o6 Ni, and energy deposition from a Magnetar. 


4.3.1 Radioactive 56 Ni model 


As can been seen in Fig. [5] the decline rate of DES13S2cmm 
is approximately that expected from the decay of 56 Co (dashed 
line). This is suggestive of the SN energy source being the pro¬ 
duction and subsequent decay of large quantities of 56 Ni produced 
in the explosion. To test this theory we fit the bolometric light 
curve of DES13S2cmm to a model of 56 Ni, following the prescrip¬ 
tion of lArnetJ (1982), whi ch is the approximation of an homol¬ 
ogous ly expanding ejecta dChatzopoulos et aPl2009l : llnserra et alj 
I2013T) where the luminosity at any given time (in erg s -1 ) is given 
by 


L(t) = Mm e (t/Tm)2 e Ni [ 
Jo 


2 u 
ZT 

0 ' m 


,(u/T m ) 2 -u/r Ni 


du 


+ eco 


(u/T m y ( 0 -u/Tf& _ g-uAto 


du 


S-y, ( 1 ) 


where u is the (time) integration variable. The energy production 
rate and decay rate for jb Ni are em = 3.9 x 10 10 ergs -1 g -1 and 
TNi = 8.8 days, while for 56 Co these values are ec 0 = 6.8 x 10 9 
ergs _1 g _1 and tco = 111.3 days. The amount of 56 Ni produced 
in the explosion is Mm- Parameterizing the rise -time of the light 
curve, Tm (from Eqns 18 and 22 of lAmettll 19821 ) is the geometric 
mean of the diffusion and expansion timescales, and is given by 
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Here E is the explosion energy, M c j is the total amount of ejected 


constant throughout, as in 

Inserra et all ( 2013)), and B is a constant 

with value ~ 13.7 (Arnett 

r«i 

00 

On 


We use d 7 to denote the gamma-ray deposition function: the 
efficiency with which gamma-r ays a r e trapp ed within the SN ejecta. 
For this function we follow 


_ __ which is also used by 

InserraetalJ (2013j) and uses the deposition function defined in 


Colgate et alJ l 1980i) . 


S y = G[1 + 2G(1 - G)(l - 0.75G)], 


(3) 
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Table 2. Best-fitting parameters for a variety of 56 Ni model fits to the bolo- 
metric light curve for DES13S2cmm. /M.max is the maximum fraction of 
*®Ni allowed in the ejecta for each fit; for DES13S2cmm the best-fitting 
model always maximizes this value. Ejecta and 56 Ni masses are given in 
units of solar masses. There are 22 degrees of freedom for each fit. 


/Ni,max 

E(10 51 erg) 

M ej 

Mm 

t 0 (MJD) 

x 2 

0.30 

31.88 

14.77 

4.43 

56508.8 

81.5 

0.40 

14.83 

10.59 

4.23 

56510.4 

76.0 

0.50 

8.22 

8.21 

4.10 

56511.5 

71.5 

0.60 

5.08 

6.68 

4.01 

56512.4 

67.9 

0.70 

3.38 

5.62 

3.93 

56513.1 

65.0 

0.80 

2.38 

4.84 

3.88 

56513.7 

62.7 

0.90 

1.74 

4.25 

3.83 

56514.2 

60.9 

1.00 

1.32 

3.79 

3.79 

56514.6 

59.4 


where G = r 7 /(r 7 + 1.6), with the ‘optical depth’ for gamma-rays 
approximately given by 



, /O.lWT 2 X 

/ 5.53 x 10 10 \ 

T 7 ~ 

^ K ' ' 4t^j / 

Vu ej (0.1 + t/r Ni ) 2 ) 


Here v e j = ^lOEk/SAEj and i s in units of cm s -1 . We no te that 
the deposition function used in Ichatzopoulos et al 1 120091 ) has a 
functional form that is similar to that used here, but in their approx¬ 
imation t 7 only deviates from « 1 at much later epochs, resulting 
in a light-curve decay rate mirroring 56 Co and largely insensitive 
to changes in the explosion energy or ejecta mass. 

We thus have four parameters in our model: the explosion 
epoch (to), the energy of the explosion (E), the total ejected mass 
(M e j), and the amount of ejected mass which is 56 Ni (Mm). We 
define the fraction of the ejecta mass that is in 56 Ni as /ni = 
Mm /Mej. We fit our model to the data with a variety of upper limits 
on /m , varying from 0.3 to 1.0. However, we follow liiiserra et all 
d2013l) in assuming a physically motivated upper limit on /m of 0.5, 
which they base on the prevalence of intermediate-mass elements 
in their spectra of lower -2 SLSN e, and the calculations of 5 6 Ni pro¬ 
duction in core-collapse SNe bv lUmeda & Nomotol ( l2008h . We re¬ 
port the best-fitting parameters and goodness of fit for these models 
in Table [2] as well as show the range of /ni models in Figure 6 to 
illustrate the scale of uncertainties in the theoretical modelling. 

We find that the 56 Ni model is not a good fit to the bolometric 
light curve, as the best-fitting model has a x 2 /dof = 2.7 (\ 2 of 59 
for 22 degrees of freedom), and is physically unrealistic (/m = 1 ). 
As can be seen in Table[2] the model prediction for Mm is relatively 
robust (3.8 — 4.4Mq), as this parameter is primarily constrained by 
the peak luminosity of the light curve. However, to match the late 
time (t > 30 days) flattening of the light curve a large amount 
of ejecta mass or low explosion energy (see Eqn. [2} is required 
to increase the diffusion time-scale, which simultaneously makes 
the fit to the post-peak decline poor while predicting a much longer 
rise time than is seen. We also note that these parameters are poorly 
constrained individually as E and M e j are correlated (Equation^, 
and this degeneracy ca nnot b e broken without the presence of high- 
quality spectral data dMazzali et alj|201~3h . 

We do not provide uncertainties on the best fit parameters for 
each /m model in Table [2] This is because the formal uncertain¬ 
ties on parameters (such as Mm) are smaller for each fit than they 
are between fits with different constraints on /m- This is due to the 
simplicity of the 56 Ni model; we do not consider herein Ni mixing, 
non-standard density profiles, nor other possible variations in this 
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Figure 6. Upper panel: Bolometric light curve of DES13S2cmm, with the 
best-fit magnetar model (green solid line) and three 5 ®Ni models with dif¬ 
ferent /Ni.max (blue lines) overplotted. The magnetar model is a better 
match to the data than the best-tit 56 Ni, and has a significantly better x 2 
than the models with a more realistic /Ni.max (see text). However both 
models have difficulty reproducing the peak luminosity, the post-peak de¬ 
cline, and the late-time flattening. Middle and Bottom panels: Temperature 
and radius evolution of the best-fitting blackbody to each photometric epoch 
which is used to construct the bolometric light curve. 


basic model that could yield different light-curve shapes. For exam¬ 
ple, it has been noted that a two-component model that is denser in 
the inne r com ponent would produce a quicker flattening post-peak 
(Maeda et al] |2003l ). 


4.3.2 Magnetar model 


As an alternative to the radioactive decay of 56 Ni, we also fit 
our derived bolometric light curve for DES13S2cmm to a magne¬ 
tar model, which proposes that SLSNe are powered by the rapid 
spinning-down of a neutron star with a n extreme magnetic field 
dWooslevll20ld : iKasen & Bildstenll201(]h . In this model the input 
energy is defined by the initial spin period (P ms , in milliseconds) 
and magnetic field strength (P14, in units of 10 14 Gauss) of the 
magnetar, and the rise-time parameter (r m , Eqn. [2} reflects the 
combined effect of the explosion energy (E), the ejected mass 
(M e j), and the opacity (+). We use the semi-analytical model out¬ 
lined in llnserra et al.l fl20 1 3T) . The luminosity of the magnetar model 
at time t is given by (in erg s _1 ) 


L(t) = 4.9 x 10 46 e“ (t/Tm)2 <5 


f 


2u 


,( U /Tm ) 2 


Bn P, 


-4 

ms 


(1 + u/r p )' 


■du, 

(5) 


where u is the (time) integration variable, t p is the spin-down 
timescale, t p = 4.7 P^ a days, and S is the deposition func- 
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tion for the magn etar radiation ; here we assume full trapping (i.e., 
<5=1) following llnserra et all 120131) . We do not include any addi¬ 
tional contribution from j6 Ni decay in this model. 

In Fig. [6] we include the best-fitting magnetar model to the 
bolometric light curve of DES13S2cmm. The best-fitting parame¬ 
ters of this model are a magnetic field of 1.43 x 10 14 Gauss, an ini¬ 
tial spin period of 5.28 ms, a diffusion timescale of 22.9 days, and 
an explosion date of MJD 56511.0. This model has a x 2 /dof = 2.0 
(X 2 of 44 for 22 degrees of freedom) and is therefore a significantly 
better fit to our data than the 56 Ni model. However, it is still clear 
that the model does not fit well : it is unable to reproduce the factor 
of two drop in luminosity post-peak within 30 days of explosion 
(although it does a better job than o6 Ni). Intriguingly the late-time 
epochs seem to fit the magnetar model well. Our time coverage 
of DES13S2cmm allows us to determine whether either model is 
able to reproduce detailed features of the light curve, but discern¬ 
ing between the two models would be improved by an even longer 
time-scale over which the two models cannot mimic one another; 
DES alone will unlikely provide longer light curve coverage as its 
observing season is typically not much more than 5 months a year. 


5 CONCLUSIONS 


In this paper we have presented DES13S2cmm, the first 
spectroscopically-confirmed superluminous supernova (SLSN) 
from the Dark Energy Survey. Using spectroscopic data obtained 
from the ESO Very Large Telescope, we measured a redshift of 
2 = 0.663 ± 0.001, and assigned a classification of SLSN-I. How¬ 
ever, we cannot exclude the possibility that DES13S2cmm is a type 
R (Fig.[4]and[5} if this is in fact a distinct SLSN type. 

Using this redshift and correcting for Milky Way extinc¬ 
tion, we find the rest-frame [/-band absolute magnitude of 
DES13S2cmm at peak to be M? r eak = — 21.05!q'q 9, consistent 
with the SLSNe definition of lGal-Yaml d2012h . Like other SLSNe 
dLunnan et al.l20l3) . DES 13S2cmm is located in a faint low stellar- 
mass host galaxy (log(M/MQ) = 8.9 ± 0.3) with sub-solar metal- 
licity. 

In Fig. we compare the bolometric light curve of 
DES13S 12cmm to fourteen similarly well-observed SLSNe-I light 
curves in literature, and see that DES 13S2cmm has one of the slow¬ 
est declining tails (beyond +30 rest frame past peak) of all SLSNe- 
I studied herein, as well as likely being the faintest at peak. We 
further find that the dispersion between the bolometric light curves 
of all SLSNe-I shown in Fig. [5] has a minimum of 0.29 magni¬ 
tudes around +30 days past peak (rest frame). This reduces to 
0.20 magnitudes around +25 days if we remove the four SLSNe- 
I above 2 ~ 1 from the measurement because of the increased 
uncertainty in their rest-frame UV luminosities. This observation 
raises t he tantalising possibility of ‘standardising’ the se SLSNe like 
SNe la dOuimbv et alj|201 llllnserra & Smartd l2014l) and starting a 
new era of supernova cosmology (King et al.(2014). Further study 
is required to confirm if SLSNe can be standardised and investigate 
further the k-corrections of SLSNe, which would be needed to use 
SLSNe as robust distance indicators. 

In Fig. [6] we fit the bolometric light curve of DES13S2cmm 
with two possible models for the power source of these extreme 
ev ents; the radioactive decay of 56 Ni and a magnetar (e.g. see 
llnserraetaD 1 201 3l) . We find that the 56 Ni model is not a good 
fit to the bolometric light curve, as the best-fitting model has a 
X 2 /DoF= 2.7 with an unrealistic fraction of 56 Ni to ejected mass 
(/Ni = 1). The model does provide a relatively robust prediction 


for the overall o6 Ni mass (3.8 — 4.4M@), as this is primarily con¬ 
strained by the peak luminosity of the light curve, but is unable to 
reproduce the late time (t > 30 days) evolution of the light curve. 
The Magnetar model provides a smaller x 2 (for the same degrees 
of freedom) than 56 Ni, but is likewise not a good fit to the data as it 
is unable to fully reproduce the drop in luminosity within 30 days 
of explosion. 

In the future, DES should find many more SLSN-I. Such data 
should allow us to further test these two models, especially if we 
can measure the (rest frame) light curves to > 100 days (see 
Fig. [6] and Inserra et al. 2013). Therefore, we have begun a new 
project entitled ‘SUrvey with Decam for Superluminous Super¬ 
novae’ (SUDSS), which will enable extended monitoring of some 
DES SN fields (in addition to several non-DES fields) to over 8 
months, greatly enhancing both our ability to detect many more 
SLSNe and measure their long light curves. 
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